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ABSTRACT 



We report the first angular diameter measurements of two classical Cepheids, FF Aql and T Vul, that we have obtained with the 
FLUOR instrument installed at the CHARA interferometric array. We obtain average limb-darkened angular diameters of Old = 
0.878 + 0.013 mas and 9 L0 = 0.629 ± 0.013 mas, respectively for FF Aql and T Vul. Combining these angular diameters with the 
HST-FGS trigonometric parallaxes leads to linear radii R = 33.6 ± 2.2 R Q and R = 35.6 ± 4.4 R Q , respectively. The comparison with 
empirical and theoretical Period-Radius relations leads to the conclusion that these Cepheids are pulsating in their fundamental mode. 
The knowledge of the pulsation mode is of prime importance to calibrate the Period-Luminosity relation with a uniform sample of 
fundamental mode Cepheids. 

Key words. Techniques: interferometry, high angular resolution ; Stars: variables: Cepheids, distances, oscillations 



£h ' 1. Introduction 



The direct detection of angular diameter variations of a pul- 
sating star using optical interferometers has been achieved for 



2004 



The 



ma ny stars no w (see e.g. Lane et al. 2002; Kervella et al. 
IMerand et ai1 l2006b: D avis et al.1 12009: Lac our et all 2009). 
combination of angular diameters with radial velocity measure- 
ments allows us to accurately estimate their distances in a quasi- 
geometrical way. This independent distance determination is es- 
sential to calibrate the Period-Luminosity (P-L) and Period- 
Radius (P-R) relations of Cepheids. Direct radius measurements 
also enable us to compare the theoret ical (e.g. Neilson et all 
1201 Ol) and indirect radius estimates (e.g. lGroenewegenl2007l) . 

Having a direct determination of the diameter, we can also 
settle whether a Cepheid belongs to the fundamental mode group 
or not, because its average linear diameter should depart from 
the classical P-R relation. As different modes yield different re- 
lations, it is essential to know the pulsation mode in order to cali- 
brate the P-L relation from a homogeneous sample of Cepheids. 

We present here the first interferometric observations of the 
Cepheids FF Aql and T Vul. The former has a small ampli- 
tude (less than 0.5 mag in V), a sinusoidal-like light curve, 
and is therefore generally designated as a s-Cepheid. The s- 
Cepheids are suspec ted to be first overtone pulsators (see e.g. 
lBohm-Vitensd[l988l) . T Vul has the same pulsation period that 
FF Aql but is expected to pulsate in the fundamental mode. 

These Cepheids also belong to binary syst e ms w ith a hot- 
ter companion. From IUE spectra, lEvans et all (1 19901) detected 
around FF Aql a companion with a spectral type between A9 V 
and F3 V. This corresponds to a magnitude difference with re- 
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spect to the Cepheid AK ~ 6 . 1 -6 .5 mag. lEvans! d 1992b found that 
the companion is a A0.8 V star, also estimated from IUE spectra, 
and corresponds to a magnitude difference AK ~ 5.2 mag. 

In this paper, we focus on determining the distance and 
radius of these Cepheids using the interferometric Baade- 
Wesselink method (IBWM). We also compare the linear radii 
with those in the literature and with published P-R relations in 
order to reveal the pulsation mode. 

2. Interferometric observations 

The observations were performed in June 2010 and A ugust 201 1 
at the CHARA Array (Iten Brummelaar et al~ll2005l) in the near 
infrared K' band (1.9 < A < 2.3 /mi) wit h the Fiber Linked Unit 
for Optical Recombination (FLUOR; ICoude du Foresto et al.l 
2003; M erand et al.l l2006ab . For these observations, we used 
baselines longer than 200 m in order to resolved sufficiently well 
the two targets. The journal of the observations is presented in 
Table [TJ 

The squared visibility of the fringes has been estimated 
from the raw data using the FLUOR data reductio n software 
(ICoude du Foresto et alj|1997t IMerand et all l2006al) . based on 
the integration of the fringes power spectrum. The raw squared 
visibilities have then been calibrated using resolved calibrator 
stars, chosen from the catalogue of M erand et al.l (2005a), using 
criteria defined to minimize the calibration bias and maximize 
the signal-to-noise ratio. The calibrators used for these obser- 
vations are listed in Table [2] They were observed immediately 
before and/or after the Cepheid in order to monitor the inter- 
ferometric transfer function of the instrument. The error intro- 
duced by the uncertainty on each calibrator's estimated angular 
diameter has been propagated using the formalism developed by 
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Table 1. Calibrated squared visibility measurements of FF Aql. 



Star 



MJD 



B 
(m) 



PA 

o 



V 1 : 



FF Aql 55352.404 


0.36 


318.98 


25.1 


0.2940 ± 0.0205 


55352.433 


0.36 


312.70 


19.3 


0.3958 ± 0.0244 


55353.431 


0.58 


273.78 


-44.0 


0.4774 ±0.0151 


55355.436 


0.03 


310.34 


16.7 


0.3584 ±0.0126 


55366.374 


0.48 


218.94 


62.9 


0.6192 ±0.0175 


55779.257 


0.83 


314.15 


20.7 


0.3839 ± 0.0095 


55781.254 


0.27 


309.83 


72.1 


0.3535 ±0.0120 


55782.317 


0.51 


230.92 


0.1 


0.5813 ±0.0159 


55784.287 


0.95 


278.37 


-48.7 


0.4678 ± 0.0200 


55786.266 


0.40 


308.46 


14.4 


0.3506 ± 0.0136 


55786.290 


0.40 


304.97 


8.7 


0.3807 ± 0.0144 


55790.230 


0.28 


249.54 


76.0 


0.5141 ±0.0126 


T Vul 55779.439 


0.25 


322.81 


-4.4 


0.5905 ± 0.0227 


55779.461 


0.25 


323.82 


-9.6 


0.6059 ±0.0210 


55779.487 


0.25 


325.80 


-15.7 


0.6148 ± 0.0222 


55780.445 


0.48 


323.12 


-6.4 


0.5279 ± 0.0294 


55780.496 


0.48 


326.81 


18.2 


0.5335 ± 0.0257 


55781.403 


0.69 


283.61 


55.9 


0.6252 ± 0.0172 


55781.423 


0.69 


270.23 


50.3 


0.7054 ±0.0213 


55782.479 


0.93 


246.61 


-18.3 


0.7524 ±0.0211 


55784.382 


0.36 


274.57 


-51.5 


0.6772 ± 0.0332 


55787.478 


0.06 


326.91 


-18.5 


0.6241 ± 0.0256 



Notes. MJD is the date of the observations (Modified Julian Date), <p the 
phase of pulsation, B the telescope projected separation, PA the baseline 
projection angle, and V 2 the squared visibility. 

Table 2. Calibrators used for the observations. 



Calibrator 


my 


m K 


Spec. Type 


(mas) 


r 
O 








FFAql 






HD 172169 


6.7 


4.1 


K4III 


0.942 ± 0.012 


13.1 


HD 169113 


7.1 


3.9 


K1III 


0.816 ±0.011 


13.3 


HD 187193 


6.0 


3.9 


K0II-III 


0.893 ± 0.012 


14.0 









TVul 






HD 198330 


7.3 


3.8 


K4III 


0.913 ± 0.012 


2.6 


HD 201051 


6.1 


4.0 


K0II-III 


0.827 ±0.011 


3.6 


HD 192944 


5.3 


3.0 


G8III 


1.085 ± 0.015 


8.5 


HD 187193 


6.0 


3.9 


K0II-III 


0.893 ± 0.012 


14.5 



Notes. ni\t,m K : magnitudes in V and K bands; 6 VD : uniform disk angu- 
lar diameter in K band; y: angular distance to the Cepheid. 



(|Perrinl2003l) . The final uncertainty of the calibrated squared vis- 
ibility includes statistical (from the dispersion of the raw squared 
visibilities obtained during the observations) and systematic er- 
rors (from the error bars on the calibrators). 

For each night we observed a Cepheid, we determined a uni- 
form disk diameter based on one or several squared visibility 
measurements (Table [TJ- Each night was assigned a unique pul- 
sation phase established us ing the average date of observation 
and the lSamus et al] d2009l) ephemeris: P = 4.470916 days and 
T = 2441576.4280 for FF Aql, and P = 4.435462 days and 
T = 2441705.1210 for TVul. 



3. The interferometric Baade-Wesselink method 

The IBWM makes use of radial velocity and angular diameter 
measurements of the star during its pulsation cycle. This in- 
volves a good phase coverage and the highest possible angular 



resolution. The integral over time of the radial velocity is then 
directly linked to the angular diameter variation. 

The presence of a companion can falsify the results if its 
magnitude is close to the one of the Cepheid, and if the orbital 
effect is detectable in the spectroscopic measurements. In our 
case, these main sequence companions are faint compared to the 
Cepheids, and do not affect the photometric signals. Only the 
radial velocities are affected by the presence of the companions 
and have to be corrected. 

3. 1. Radial velocity integrations 

Among pub l ished radial velocities data of FF Aql, we chose 
lEvans et all (Il990t) for two reasons. First, the orbital motion 
causes by the presence of the companion was already re- 
moved from measurements, only leaving the displacement due 
to the pulsating photosphere. Secondly, the radial velocities 
were extracted using the cross-correlation method. As shown by 
iNardetto et al] (l2004l) . the method used can affect the distance 
determination via the choice of the projection factor (hereafter 
p-factor), which we w ill introduce in Sec | 3 .2l For T Vul, we use 
the radial velocities of iBersier et al] (ll994l) . also from the cross- 
correlation method. The amplitude of the orbital motion was not 
detected with this data so we did not apply any corrections due 
to the presence of the companion. For both Cepheids, we re- 
determined the pulsation phases using the ephemeris presented 
previously. 

The radial velocities, which were acquired at irregular 
phases, must be numerically integrated to get the radial displace- 
ment. The best way to achieve a robust integration and avoid 
noisy data is to interpolate the data before integration. For this 
purpose, the radial velocities were smoothly interpolated using 
a periodic cubic spline functi on, define d by floating nodes. This 
method was already used by Mer and et al.l (I2007I) and we refer 
the reader to this paper for a detailed explanation. The interpo- 
lated curves for our two Cepheids are presented in Fig. [TJ The 
dispersion of the residuals are cr — 61 m s _1 and <x = 68 m s , 
respectively for T Vul and FF Aql. 

3.2. Distance determination 

Our angular diameter measurements were then fitted to the radial 
displacements to get the distance d: 



0ud(T) - 0ud(O) 






(t) - Vy)dt, 



(1) 



where p is the so-called /^-factor, defined as the ratio between 
the pulsation velocity and the radial velocity (measured by spec- 
troscopy), #ud me interferometric angular diameter for a uni- 
form disk model, k the ratio between the uniform and the limb- 
darkened stellar diameter, and v y the systemic velocity. 

The parameters of the fit are the average angular diameter 
#ud(0), the quantity d/lcp and a possible variation in the pulsa- 
tion period between interferometric and spectroscopic observa- 
tions. 

The /^-factor depends on detailed modelling of stellar atmo- 
sphere, but such models are only available for a few Cepheids 
and we conventional ly choose a multiplicative constant factor. 
INardetto et al.l (120041) showed that the choice of a constant factor 
instead of a time-dependent ones gives a systematic error in the 
final distance of ~ 0.2 %, which is below the best relative pre- 
cision of current distance estimates. iMerand et al.l (l2005bl) pre- 
sented the unique measured value p - 1.27 + 0.06 for 6 Cep. 
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Fig. 1. Radial velocity of FF Aql and T Vul. The solid line is the periodic spline function defined with four adjustable floating nodes 
(large black dots). The lower panel displays the residuals of the fit. 
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Fig. 2. Uniform disk angular diameter variation of FF Aql and T Vul. The solid line is the integration of the radial velocity curve 
with the distance, average angular diameter, and period change as fitted parameters. The dashed line includes the period change. 
The lower panel displays the residuals of the fit. 



Although there is no agreement in the literature on the optimum 
value of the p-factor, it is clear that it depends on each Cepheid. 
Different authors use either constant p-factor values (ranging 
from 1.27 to 1.5), or a dependence with the pulsation period. 
The period dependence is likely to reflect the dependence with 
the limb darkening. For this work, we chos e the latest p-facto r 
relation p = 1.550*0.04 - 0.186 ±0 .o6logP of lStormet all(l201lb . 
because this is an empirical relation based on the largest sample 
of Cepheids, and covering the widest range of periods. We have 
for both stars p = 1 .43 + 0.06. 

The k parameter is generally determined assuming a limb- 
darkened disk model and is expe cted to vary wi t h the p ulsation 
phase. However, in the near-IR, iMarengo et alj (120031) showed 
from hydrodynamic models in spherical geometry that k is con- 
stant wi th the pulsation at a level of 0.2 %. From hydrostatic 
models, iMerandl d2008l) found a variation of k of the order of 
0.4 %. Our interferometric measurements are not sensitive to this 



variation and we therefore assumed a constant limb-darkening 
coeffic ient. The hyd rostatic model computed from Claret coeffi- 
cients (IClareJl2000l) . with the average stellar param eters T e g — 
6000 K , logg = 2, Fe/ H = 0.0 and V t = 4 km s~ ' from lLuck et"ai1 
d2008h for FF Aql and lAndrievskv et all d2005b for T Vul, corre- 
sponds to k = 0.985. The IBWM fit yields for FF Aql a distance 
d - 363 + 73 pc and an average angular uniform disk diameter 
6» UD = 0.865 + 0.013 mas, with a total reduced^ 2 of 0.4. For 
T Vul, we obtain d = 415 + lOOpc and 6» UD = 0.620 ± 0.0 12 mas, 
with a total reduced^- 2 of 2.5. These results are summarized in 
Tableland the angular diameter variations are plotted in Fig. [2] 
We also allowed a change in the period for FF Aql and T Vul of 
~ 14 s and ~ -26 s, respectively. 
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Table 3. Cepheid average angular diameters and distances determined through the application of the interferometric B W method. 



Star 



(mas) 



(mas) 



d 

(pc) 



"HST 

(pc) 



R 

(Re) 



FFAql 0.865 ±0.013 0.878 ±0.013 363 ± 73 356 ± 23 33.6 ± 2.2 0.4 
TVul 0.620 ±0.012 0.629 ±0.013 415 ± 100 526 ± 64 35.6 ±4.4 2.5 

Notes. The linear radius R were estimated using the more accurate distances rfnsT derived by Be nedict et al.l d2007h . 



3.3. Linear radius 

We transformed the uniform disk diameters to limb-darkened 
disk diameters using the value of k previously chosen, we get 
LD = 0.878+0.013mas for FFAql and LD = 0.629+0.013 mas 
for T Vul. These measured values, which have an accuracy 
~ 1.5 %, are in good agreement with those predicted from 
IR surface brightne ss method (e.g. IMoskalik & Gorynya 2005; 
lGroenewegenll2007l) . The fitted distance of FF Aql is also con- 
sistent with the derived value d = 35 6 ± 23 pc from the HST- 
FGS parallax of iBenedictet al.l d2007l) . The same authors give 
for T Vul d - 526 + 64 pc, that is ~ 27 % larger than our va lue, 
but within 1.1 cr. As the values from iBenedict et all d2007l) are 
more accurate and do not depend on a p-factor assumption (com- 
pared to our fitted distances), we preferred them to evaluate the 
linear radius. We therefore used the HST-FGS distances in the 
BW fit and obtained R = 33.6 + 2.2fl Q and R = 35.6 + 4.4 R Q , 
respectively for FF Aql and T Vul. 

In our analysis, we did not take into account a possible cir- 
cumstellar envelope (CSE) emission, that could lead to an over- 
estimate of the angular dia meter. From the spectra l energy dis- 
tribution of FF Aql given bv lGallenne et al.l (1201 ll) . the infrared 
excess caused by the CSE appears around 10 /mi, while it is neg- 
ligible at 2.2 /mi. We assumed this is also the case for T Vul. 



3.4. p-factor 

We can use the HST-FGS distances in Equ.[TJto directly fit the 
p-factor, as already done bv lMerand et al.l(l2005bl) for 6 Cep. We 
get p = 1.40 ± 0.28 and p = 1.81 ± 0.44, for FF Aql and T Vul 
respectively. These estimates can be compared wit h the empiri- 
cal an d the oretical p-P rel a tions, respectively from lStormet al.l 
d201ll) and iNardetto et all (120091) . For FF Aql, the p-factor is 
compatible with the value derived from the empirical relation, 
while the estimate of T Vul is within l<x. Although at odds with 
the theoretical relation, our estimates are not accurate enough to 
constrain it. 

A direct determination of the p-factor is a powerful con- 
straint to numerical atmosphere models of Cepheids, however 
our work does not contain enough measurements to have an ac- 
curate estimate of this parameter for these stars. 

4. Period-Radius relation 

The estimated linear radii can be compared with empirical and 
theoretical Period-Radius (P-R) relations, that are of the form 
log R = a log P + b . We ch ose the recent theoretical calculations 
from iNeilson et al.l d2010h . who used a pulsatio n-driven mass 
loss m odel. We also selected the recent works of Groenewegen 
d2007l) . who used Cepheids with known distances and angular 
diameters to derive an empirical P-R relation. These two rela- 
tions are only valid for fundamental mode Cepheids. 

In principle, the position of a Cepheid in the P-R diagram 
can reveal its pulsation mode, if the estimated radius is suffi- 



ciently accurate. We plotted in Fig. [3] the P-R relations (thick 
curves) and our estimated linear radii. The intrinsic dispersion 
of each relation is represented with thinner curves (with the 
same colour and symbol). We notice that both Cepheids are con- 
sistent (within the intrinsic dispersion) with stars pulsating in 
a fundamental mode. We also plotted with triangles the radii 
assuming an overtone pulsation, using the ratio of the over- 
tone period Pi to the fundamental period Pq given by Pi/Po = 
0.720 - 0.027 logPo (lAlcocketal.lll995l) . In addition, we rep- 
resented with squares the radii of the star assuming first over- 
to ne radii, by using the theoretical first overtone P-R relation 
of iBono et alJ (1200 ll) . The overtone mode seems not compatible 
with the P-R relations, and therefore suggest that FF Aql and 
T Vul are pulsating in the fundamental mode. 
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Fig. 3. Period-Radius relation for fundamental mode Cepheids. 
The colour curves represent empirical or theoretical P-R rela- 
tions from the literature. Each P-R intrinsic dispersion is plotted 
with a thinner curve, with the same colour and symbol. The dots 
are our estimated linear radii, the triangles are converted for first 
overtone pulsation, and the squares represent the theoretical first 
overtone radii. 



5. Conclusions and discussion 

We report the first interferometric observations of two classi- 
cal Cepheids, FF Aql and T Vul, with the CHARA/FLUOR in- 
strument. The angular diameter variations were monitored over 
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the pulsation cycle and combined with previously published ra- 
dial velocity measurements to get an independent estimate of the 
distance and mean angular diameter. The precision achieved is 
~ 2 % on their average angular diameter estimates. 

We found for FF Aql a good agreeme nt (0.3cr) between ou r 
fitted distance and the HST-FGS estimate dBenedict et all2007l) . 
For T Vul, a discrepancy is found between the two estimates, our 
value being ~ 27 % smaller, but within 1.1 cr. The main source 
of bias in the use of the BW method is the choice of the /^-factor, 
and that could explain the discrepancy encountered. We can also 
note that more interferometric measurements, with higher ac- 
curacy, would also reduce this discrepancy. To avoid a possible 
bias on the distance due to the /^-factor assumption, we coupled 
our angular diameters with existing trigonometric parallaxes to 
evaluate the linear radii. 

When we compared the linear radii with empirical and the- 
oretical P-R relations, we showed that these two short period 
Cepheids are pulsating in their fundamental mode. Our work 
confirms the pulsation mode of T Vul, alr eady flagged as fun- 
damental pulsator by several authors (e.g. iGroenewegenl 12000, 
and references therein). However this is not the case for FF Aql. 
This s-Cepheid, with small light amplitude and practically sinu- 
soidal light c urve, is usually considered as pulsating in the first 
overtone Ce.g. lSachkovll997l) . As the pulsation mode is generally 
identified by Fourier decomposition of the light curve or the ra- 
dial velocity curve, the presence of a companion around FF Aql 
could lead to an error in the determination of the pulsation mode. 
If the difference between the magnitude of the Cepheid and that 
of its companion is small, the apparent magnitude of the Cepheid 
will be falsified. The radial velocity measurements can also be 
altered because of the orbital effect. 

The knowledge of the pulsation mode is of prime importance 
to calibrate the P-L relation with an uniform sample of funda- 
mental mode Cepheids, since different modes will yield different 
relations. 
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